Heterogeneous condensation of the Lennard-Jones vapor onto a nanoscale seed particle 
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The heterogeneous condensation of a Lennard-Jones vapor onto a nanoscale seed particle is studied 
using molecular dynamics simulations. Measuring the nucleation rate and the height of the free 
energy barrier using the mean first passage time method shows that the presence of a weakly 
interacting seed has little effect on the work of forming very small cluster embryos but accelerates 
the rate by lowering the barrier for larger clusters. We suggest that this results from a competition 
between the energetic and entropic features of cluster formation in the bulk and at the heterogeneity. 
As the interaction is increased, the free energy of formation is reduced for all cluster sizes. We also 
develop a simple phenomenological model of film formation on a small seed that captures the general 
features of the nucleation process for small heterogeneities. A comparison of our simulation results 
with the model shows that heterogeneous classical nucleation theory provides a good estimate of 
the critical size of the film but significantly over-estimates the size of the barrier. 
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I. INTRODUCTION 

The heterogeneous condensation of a vapor onto a sub- 
strate is a key step in a wide range of chemical and 
physical processes that occur in both nature and technol- 
ogy. For example, dust and pollutant aerosol particles, 
ranging in size from several microns down to just a few 
nanometers, serve as cloud condensation nuclei in the at- 
mosphereii^, while nanoscale structured surfaces provide 
templates for the controlled nucleation and growth of 
complex materials such as colloidal metamaterials-. The 
general principles of heterogeneous nucleation on bulk 
surfaces are well knowr^, but the details of how the 
process changes as the size of the heterogeneity becomes 
microscopic, and is the same size, or smaller, than the 
nanometer-sized fluctuations involved in the nucleation 
process itself, are not well understood. Recent experi- 
ments^ have shown that vapor condensation onto small 
clusters and nanoparticles, with diameters of 1 — 24nm, 
occurs at supersaturations well below those predicted by 
classical nucleation theory^"— (CNT) and the Kelvin re- 
lation, which suggests very small particles are far better 
at activating nucleation than expected. Similarly, Sear- 
found that a microscopic heterogeneity, consisting of a 
single fixed spin, was sufficient to accelerate the nucle- 
ation rate in a two-dimensional Ising model by as much 
as four orders of magnitude. This raises interesting ques- 
tions regarding how small a heterogeneity can be and 
still have an influence on nucleation and whether there 
are practical limits to our ability to study homogeneous 
nucleation in the presence of small concentrations of im- 
purities. 

The work of forming a critical size cluster is the central 
ingredient of CNT. The phenomenological expression for 
the work of forming an n-sized embryo of the new phase 
can be generally written as 

2/3 



AG'(n) = nA^ 



(1) 



where A/i is the difference in chemical potential between 
the bulk stable and metastable phases and a is a geo- 
metric constant that is dependent on the shape of the 



embryo and proportional to the interfacial surface ten- 
sion, 7. The first term represents the free energy gained 
by moving particles from the metastable mother phase to 
the more stable phase and is negative, while the positive, 
second term accounts for the free energy cost of introduc- 
ing an interface between the two phases. The functional 
form of Eq. [T] highlights the activated nature of the nu- 
cleation process and shows that embryos must overcome 
a free energy barrier, AG(n*), associated with making 
a critical cluster of size n* , before they can grow spon- 
taneously into the new phase. CNT then expresses the 
rate of nucleation as J = Aexp[AG(n*)/fcr], where the 
critical barrier height is obtained from dG{n)/dn = 0, A 
is a kinetic prefactor that depends on the critical clus- 
ter size and the curvature at the top of the barrier, k is 
Boltzmann's constant and T is the temperature. 

The simplest phenomenological model used in CNT 
to describe the condensation of a supersaturated vapor, 
which involves the capillarity assumptions of uniform 
bulk densities for both phases and a sharp interface char- 
acterized by the bulk, planar surface tension, gives the 
critical barrier and critical cluster size for the homoge- 
neous nucleation as 
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respectively. Here, vi is the volume per molecule in the 
bulk liquid phase, S is the supersaturation and we have 
used A/i sa kT\x\S as well as assuming that the nucleus 
is spherical in shape. 

In the presence of a macroscopic surface, nucleation 
can occur via a heterogeneous pathway. The nucleus 
then forms a droplet, contacting the wall with an an- 
gle, 9, that characterizes the interaction of the material 
with the wall and is related to the wall-liquid (7^1/) and 
wall- vapor (7^,,,) surface tensions through Young's equa- 
tion, 7t,i cos 6 = 7^1, — 7u,/ . The extension of CNT to 
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heterogeneous nucleation by TurnballiS. gives 

AGhetin*) ^ fW^Gi^omin*), (4) 

where < f{9) < 1 is solely a function of the contact 
angle. Eq|4] shows that the barrier is always reduced be- 
cause the wall contributes a portion of the interfacial free 
energy. In the limit where the contact angle goes to tt, the 
liquid is no n- wetting and f{9) tends to one. Nucleation 
then occurs homogeneously in the bulk. If the liquid and 
surface are highly attractive so the liquid completely wets 
the wall, 9^0 and f{9) — >■ 0, causing the barrier to go to 
zero. Despite the extremely simplified approach involved 
in CNT, recent simulation studies of heterogeneous nucle- 
ation in hard spheres colloidsi2^J^5. and the ising modell^ 
have shown that Eq. H] is generally correct, as long as 
an additional term involving the three-phase contact line 
tension is included in AG{n). CNT also suggests that as 
the heterogeneity becomes more microscopic its ability 
to activate nucleation is greatly reducedii. 

The goal of the present work is to quantify the effect 
of a very small seed particle on the condensation of the 
Lennard - Jones (LJ) vapor, by measuring both the nucle- 
ation rate and free energy barrier associated with hetero- 
geneous nucleation, and compare the results with a sim- 
ple thermodynamic model. Homogeneous nucleation in 
the LJ system has been studied by simulationi^""— as 
has heterogeneous nucleation on macroscopic surfaces^^. 
However, studies involving small, seed particles in the 
LJ vapor have focused on examining the qualitative ef- 
fects of highly attractive seeds, i.e. where the vapor - 
seed interaction is in the order of being ten times more 
attractive than the vapor-vapor interactioiiSi. By study- 
ing the impact of seeds with interactions very similar to 
those of the vapor phase, we hope to understand how 
small heterogeneities perturb the bulk system. We can 
then directly address the question as to whether small 
impurities can seriously impact our ability to measure 
homogeneous nucleation rates, as well as assess the abil- 
ity of CNT to predict heterogeneous nucleation on seed 
particles. The paper is organized as follows: Section 2 
describes the development of a simple thermodynamic 
model that provides the appropriate CNT comparison 
for our nucleation rates and barriers calculated, using 
molecular dynamics (MD) simulations, in Section 3. Our 
discussion is contained in Section 4. 



the particle giving rise to a film-particle composite of ra- 
dius r2 and leaving rii = N—n2 atoms in the vapor phase 
which is treated as an ideal gas. At constant N, V, T, the 
Helmholtz free energy, F, is the appropriate thermody- 
namic potential and variations in F are given by 

dF = dU - TdS , (5) 

where, U — Ui + U2, is the total internal energy and 
S — Si + S2 is the total entropy. Here, the subscripts 1 
and 2 represent the vapor and liquid phases respectively. 
Variations in the U are given by 

dUi — TdS — pidVi + /iidrii , (6) 

and 

dU2 = TdS2 - P2dV2 + fJ,2dn2 + 712^^12 -I- 702^^02 , (7) 

where pi and p2 are the respective pressures of each 
phase. We will assume that the interfaces between the 
phases are sharp, consistent with the capillarity approx- 
imation, so that the area of the liquid-vapor interface 
is given by, A12 = 47rr|, and the solid-liquid surface 
area, Aq2 — Atttq. 712 and 702 are the planar sur- 
face tensions of the interfaces respectively. Substitut- 
ing Eqs. [6] and [7] into Eq. [5] while using the conserva- 
tion laws, dVi — —dV2 and drii = —dn2, along with 
dAi2 = 2dV2/r2 and dAoi = yields, 

dF = ~{p2 ~pi- 2ji2/r2)dV2 + (/i2 - Aii)d?i2 ■ (8) 




II. THERMODYNAMIC MODEL 

We begin examining the problem of heterogeneous nu- 
cleation onto a nanoscale seed particle by considering the 
simple thermodynamic model described in Fig. [1] which 
represents the heterogeneous extension of the liquid drop 
mode l^^i^^ . The system has a fixed total number of par- 
ticles N, a fixed volume V, and a constant temperature, 
T. The heterogeneity is considered to be an insoluble, 
spherical, solid particle with radius rp, while the liquid 
forms a uniform film of n2 atoms that completely wets 



FIG. 1. N,V,T model consisting of an insoluble heterogeneity 
of radius ro completely wet by a thin film of 722 liquid atoms 
giving the film-seed composite a total radius of r2. The vapor 
is an ideal gas of ni atoms. 

At equilibrium, dF = 0, and Eq.[S] yields the expected 
equality of chemical potentials between the two phases 
and the Laplace relation for the difference in pressures 
across a curved interface, 

M2 = Ml > (9) 
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Ap = p2-pi = 27i2/r2 



(10) 



A more detailed description of coexistence can be ob- 
tained by using the Gibbs-Duhem relation for each phase. 
For the vapor, we have SidT — Vidpi + nidpLi = 0, which 
leads to 



Pi 

eq 

Pi 



(11) 



upon integration at constant T and where we have chosen 
the reference pressure, p^', as the coexistence pressure of 
the vapor in contact with a liquid, with a planar interface. 
The Gibbs-Duhem relation for the liquid and its associ- 
ated interface is, S2dT — V2C?P2 + n2dfi2 + ^12^712 = 0. 
If we assume that the surface tension is independent of 
pressure and that the liquid is incompressible, the inte- 
gration of the Gibbs-Duhem relation gives, 



M2(P2) - Ml bi ) = V2iP2-Pi ) 



(12) 



where V2 is the volume per molecule in the liquid phase. 
We have also taken advantage of Eq. [9] and the fact that 
Eq. 1101 gives us = Pi' i'^ the limit r2 00, i.e. at 
a planar interface. Combining Eqs. [TUl [11] and [T^] yields 
the coexistence equation. 
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{N - n2)kT 
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(13) 
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Eq. [12] is equivalent to the Kelvin relation and can be 
solved to find the equilibrium size of the film. In an open 
system, there would be one solution corresponding to the 
unstable equilibrium of the critical-sized film. Films that 
are thinner than the critical size tend to evaporate while 
thicker films grow spontaneously into a macroscopic drop 
as they consume molecules from a continually replenished 
vapor. However, in our closed system, N is fixed and the 
growth of the liquid film depletes the number of molecules 
in the vapor phase, causing the supersaturation to de- 
crease so that the film and vapor must eventually come 
into stable equilibrium. 

Fig [5] shows the equilibrium film sizes obtained from 
Eq. [T3]as a function of the total volume of the system for 
different sized seed particles, where we have used the LJ 
parameters of Baidakov et al^^ for Peq,V2 and 712. All 
the systems exhibit a or evaporation volume, above 
which it is not possible to stabilize a film of any size and 
the stable equilibrium state consists of a dry particle sur- 
rounded by vapor. In the limit that tq — ^ 0, we recover 
the solutions describing the formation of a liquid drop 
in the N, V, T ensemble^^ and they have been included 
here for comparison with the heterogeneous nucleation 
case. The liquid drop model has two solutions below the 
evaporation volume. The large n2 cluster belongs to the 
stable droplet in equilibrium with the vapor while the 



small droplet is the unstable, critical sized droplet that 
must be formed before the droplet can grow. We note 
two important features: firstly, there is always a bar- 
rier to the formation of the droplet, even at very small 
system volumes, which correspond to high initial super- 
saturations. Secondly, the droplet size remains finite at 
the evaporation volume. When a small insoluble particle 
is introduced into the system, we see that there are still 
two solutions at volumes just below the evaporation vol- 
ume, but now the size of the small critical cluster tends 
to zero as the volume of the system is decreased so that 
there is no barrier to film formation at small enough vol- 
umes. As the size of the solid particle is increased, we 
eventually reach a seed size for which there is never a 
nucleation barrier associated with the formation of the 
film and the stable film grows continually from ^2 = 
at the evaporation volume. This final result is consistent 
with the results of heterogeneous nucleation of films on 
bulk surfaces that exhibit complete wetting. 




FIG. 2. Coexistence solutions for the film size, 712, from 
Eq. [131 with = 512 and T* = 0.67, in the presence of 
different size heterogeneous particles, ro/a = 0,2.0 and 3.8. 
The locally stable and unstable solutions are represented by 
the solid and dashed lines respectively. 

Starting from Eq. [5] using the expressions for the 
chemical potential from Eqs. [TTl and [12] and noting that 
dV2 = V2dn2 for a incompressible liquid, gives us an equa- 
tion for dF in terms of dn2 that can be integrated at 
constant temperature to yield. 



AF{n2) - ^^(^2) - ^^(0) 
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(15) 
(16) 



where po = NkT/{V — (4/3)7rrQ) is the pressure of the 
initial vapor before any film is formed. Eq. [T^] repre- 
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sents the work of forming a film of n2 molecules starting 
from an infinitely thin film with n2 — that wets the 
particle. In principle, there is also an additional term, 
47rrQ(7i2 — 7oi), which represents the free energy change 
associated with the wetting process, but this will be a 
negative constant and is not included here as it does not 
effect the probability of observing the critical sized clus- 
ter within the context of the model. Plots of Eq. in 
Fig. [3] for different volumes confirm that the small n2 
solutions represent unstable equilibrium solutions, while 
the large n2 films are locally stable. They also highlight 
the global stability of the films and show that there is a 
range of volumes for which the thick film is locally stable, 
but metastable with respect to the dry heterogeneity. 



H 

5^ 




FIG. 3. AF{n2)/kT as a function of n2, with = 512, 
ro/cr = 1.0 and T* = 0.67, for different system volumes. 



III. MOLECULAR DYNAMICS SIMULATION 
STUDIES 

Molecular dynamics simulations in the canonical, 
(iV, y, T) ensemble are now used to study the conden- 
sation of a vapor onto a nanoscale heterogeneity. We 
model the composite system (vapor -I- heterogeneity) as 
a cut, but unshifted, Lennard - Jones mixture with the 
interaction potential. 



(17) 



where e^j and cr.y represent the energy and length inter- 
action parameters. The vapor and heterogeneity are de- 
noted as component land 2 respectively. We use N — 512 
vapor particles with en = 1.0, and a single heterogeneity 
particle with interactions in the range ei2 = 1.0 — 3.0. 
This system size has been shown to be large enough that 
the depletion of the vapor phase caused by nucleation 
in the N, V, T ensemble is small so the barriers are the 



same as those calculated in an open system21. All parti- 
cles have the same size (an = (722 = o'i2 =1), the seed 
and gas particles have the same mass, m, and the poten- 
tial is cut at Tc = 6.78(jy , which should be long enough 
to approximate the full LJ potential and is consistent 
with a previous study of homogeneous condensation^^. 
The simplicity of this model allows us to treat the het- 
erogenous seed particle as another atom in the molecular 
dynamics simulation so it is free to translate throughout 
the system. 

Simulations are carried out using the Gromacs Pack- 
age2^, with the leap-frog integration scheme. The veloc- 
ity rescaling thermostat is employed to maintain the sys- 
tem at a reduced temperature, T* — kT/en = 0.67, as 
this provides an efficient method for temperature control 
that does not appear to significant infiuence the kinetics 
of nucleation, even though the particle dynamics are per- 
turbed in a non-physical way2S. The volume of the sim- 
ulation cell is chosen to ensure the initial starting condi- 
tions correspond to a particular supersaturation, defined 
as S" = p/peq, where p is the vapor pressure and p^q is the 
equilibrium coexistence pressure at T as given by Ref.— 
Time is measure in reduced units, t* = t{eii/my^^ /an 
and we use periodic boundary conditions. 

The nucleation rate and free energy barrier are ob- 
tained using the mean first passage time (MFPT) ap- 
proach introduced by Reguera et a P°i'^° . For each state 
point studied, we obtain 200 initial starting configura- 
tions by simulating the vapor phase in the absence of 
the heterogeneity for 10^ time steps, saving configura- 
tions every 5000 time steps. The heterogeneity is in- 
serted randomly into the vapor, but ensuring that it is 
not placed within a of any vapor molecule. The MD 
trajectory is then followed as a function of time and the 
cluster size distribution is analyzed every 1000 time steps 
until the system is nucleated. Clusters are identified us- 
ing the Frenkel cluster criteriai^, which identifies liquid- 
like atoms as those particles that have at least five other 
atoms within a distance of 1.5f7, and considers two liquid- 
like atoms within a distance of 1.5a to be in the same 
liquid cluster. As a computational convieniece, we allow 
the seed particle to count as a neighbor when identifying 
liquid particles. 

The MFPT for each cluster size, r(n), is obtained by 
measuring the first time the largest cluster in the sys- 
tem reaches the size n in a simulation trajectory, and 
averaging over the 200 trajectories. However, as a result 
of the intrinsic nature of the cluster dynamics and the 
fact that we only sample configurations periodically, the 
cluster growth is non-nionotonic in time and does not 
proceed through a series of single particle additions or 
losses. It then becomes necessary to correct the time at 
which an n-sized cluster is observed for the first time in 
the simulation if a given cluster is missed in a trajectory 
all together, or if it is sampled out of order, i.e. when a 
small cluster is observed only after a larger cluster has 
already been sampled. We achieve this by assigning any 
small clusters that have not already been seen in a given 
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trajectory, the same time that is assigned to the next 
largest cluster when it first appearsS-^. 

When the barrier is high enough for the steepest de- 
scent approximations to hold, the MFPT can be repre- 
sented by^, 



TABLE I. Summary of simulation results. 



T{n) 



2 



[1 + erf{{n — n*)c)] 



(18) 



where n* is the critical size, tj is the nucleation time, 
which is related to the steady state nucleation rate as 
J = 1/tjV, c is associated with the Zeldovich factor, Z, 
as c = y/nZ and erf{x) — dx is the error 

function. This allows us to extract J, n* and c by using 
them as parameters in the fit of Eq. [TH] to our data. The 
free energy barrier is given by^i, 



PAF{n) 
where 
B{n) = - 



l3AF{no)+ln 



B{n) 
B{no) 



PstH 



Pst{n')dn' 



dn' 



-ib) - Tin) 



r{h) 



(19) 



(20) 



Here, h is the upper absorbing boundary, taken as 6 = 60 
in this work, while no = represents the lower reflecting 
boundary and reference state. Pst{n) is the steady state 
probability that the largest cluster in a given configura- 
tion from the ensemble of runs is of size n. 

A full summary of the nucleation rates, barrier heights 
and critical cluster sizes for all state points studied can 
be found in Table ID Fig. U shows T{n) for systems 
with an initial vapor pressure of S = 10.43, for differ- 
ent seed atoms with ei2 = 1.0, 1.5, 2.0 while Fig. [5]shows 
the free energy barriers. To estimate the error in our 
calculations, we divide the ensemble of 200 runs into 
10 blocks and calculate r(ri) and AF(n)/kT for each 
block. The error bars in T(n) represent the standard 
deviation of the block averages. The standard deviation 
in AF{n)/kT w 0.5 near n* , but this grows to between 
1-2 for larger cluster sizes. With ei2 — 1.0, the hetero- 
geneity is simply an additional vapor molecule and we 
obtain J = 0.04 x 10^^cm~^s~^ and n* — 20 from the 
MFPT, where we have used the LJ argon parameters, 
a = 0.3405nm, e/k = 120 K and m = 6.631 x 10"^^ 
kg to make the conversion from reduced units. The bar- 
rier calculations give AF{n)/kT = 8.78 and n* = 19. 
While our rate is approximately 33% higher, and our 
barrier approximately 0.5fcT lower, than that obtained 
by Wedekind et al22., the results are comparable when 
the slight difference in the definition of the supersatu- 
ration and the error in our calculations are taken into 
account. 

The general effect of increasing the attraction of the 
seed particle is to increase the rate by lowering the barrier 
and decreasing the critical cluster size. With ei2 = 1.5, 
the work of forming small clusters with sizes 1-5 is the 
same as in the homogeneous case, with no seed, but then 
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the heterogeneous barrier becomes lower for larger cluster 
sizes. To understand this, we analyze the MD trajecto- 
ries and identify both the largest liquid cluster, which is 
our order parameter, and the largest cluster containing 
the seed. Fig. [5] shows the two measures do not always 
coincide and most of the fluctuations involving small clus- 
ters do not contain the seed so the probability of seeing 
these clusters in the ensemble remains unchanged by the 
presence of the seed. However, the seed is always part 
of the cluster that eventually fluctuates over the barrier, 
suggesting that the added attraction is sufficient to help 
build the larger clusters and make them more probable 
that they would normally be, causing the barrier to de- 
crease. 




50 60 



FIG. 4. r(n) as a function of n, for systems with TV = 512, 
S — 10.43 and a seed particle with ei2 — 1.0 (no seed), 1.5 
and 2.0. The points represent data obtained from simulation 
and the error bars are the standard deviation of the block 
averages. The solid lines are best fits of Eq. [l8]to the data. 

Some of the fiuctuations in Fig. [5] seem to suggest 
that the largest cluster liquid containing the seed inter- 
mittently drops to zero while the largest liquid cluster 
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FIG. 5. AF{n)/kT as a function of n, for systems with A'^ = 
512, 5* = 10.43 and a seed particle with ei2 = 1.0 (no seed), 
1.5,2.0 and 2.5. 



FIG. 6. Size of the largest cluster (squares), and size of the 
largest cluster containing the seed (joined circles) as a func- 
tion of time in a MD trajectory with ei2 = 1.5 and S = 10.43 



in the system grows. Following the system using the 
Stillinger cluster criteria^, which counts particles within 
a distance of 1.5(t of each other as being part of the 
small clusters, allows us to identify those atoms that are 
loosely bound to the liquid-like particle, but may not 
have enough neighbors to be liquid-like themselves. We 
find the largest Stillinger cluster in the system and the 
largest Stillinger cluster containing the seed are identi- 
cal for the larger clusters going over the barrier, but not 
when the small clusters are fiuctuating in the metastable 
vapor. This emphasizes that the fluctuations involving 
small clusters are independent of the seed but that the 
seed is always connected to the largest liquid cluster as it 
grows, even if it is not always wet by five or more neigh- 
bors. 

The formation of monomers, dimers and trimers still 
appears to be relatively independent of the seed with 
ei2 = 2.0, (Fig.[7l), but the seed appears in most, but not 
all, of the larger fluctuations, and this is reflected in the 
barrier. For example, at S = 10.43, the seed is observed 
to be connected to trimers and 5-omer 55% and 79% of 
the time, respectively. The seed is found to be in clus- 
ters larger than n = 10 more than 96% of the time. As 
the supersaturation is decreased, the fluctuations of the 
bulk vapor naturally become rarer and most clusters are 
connected to the heterogeneity. Finally, when ei2 = 2.5, 
the nucleation times are very rapid and we are probably 
reaching the limits at which the method can be applied. 

We can further quantify the location of the seed par- 
ticle within the larger clusters by calculating ps{r), the 
radial probability of finding the the seed a distance r 
from the centre of mass of the cluster, for different n-size 
clusters, and comparing this with the density profile of 
the cluster (see Fig. |5]). While the small cluster sizes of 
10 and 14 give rise to broad density profiles where the 
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FIG. 7. Size of the largest cluster (squares), and size of the 
largest cluster containing the seed (joined circles) as a func- 
tion of time in a MD trajectory with ei2 = 2.0 and S = 10.43 



density of the core is very much lower than the bulk, it 
is clear that the seed is generally located in the center of 
the cluster, so that it is effectively wet by the condensing 
vapor. 

To study the effects of supersaturation on heteroge- 
neous nucleation, we calculate l^F{n)/kT for a system 
containing a seed particle with ei2 = 2.0 over a range of 
S (Fig. El) and compare the size of the critical nucleus 
and height of the barrier with our CNT-based model in 
Fig. [TOl The supersaturation for the model is defined 
as po/peq in order to maintain the self-consistency. This 
assumes that the LJ vapor phase is well described by 
an ideal gas at the temperatures and densities studied, 
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(l/ln(S))' 



FIG. 8. The density profile (solid line) and ps(r) (dashed line) 
for n=10 (top) and n=14 (bottom) clusters, with S — 10.43 
and ei2 = 2.0. 



which is generally true away from the critical poin t^^i^^ , 
but still leads to us underestimating the supersatura- 
tion in the model by approximately 14% compared to 
the simulation at S" = 10.43. Our results for hetero- 
geneous nucleation essentially mirror those observed in 
homogeneous nucleation, in that the critical cluster size 
predicted by CNT is close to our simulation results, but 
the nucleation barriers are overestimated by a 100% or 
more, which would lead to many orders of magnitude 
error in the rate. 
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FIG. 9. AF{n) /kT as a function of n at different supersatu- 
rations for a system with TV = 512 and a seed particle with 
£12 = 2.0 



FIG. 10. (Top) Critical cluster size, n* as a function of 
(l/lnS)'^. The circles and squares were obtained from fits of 
Eq. [T8]to our simulation values of r(n) and the AF{n) /kT 
curves in Fig [51 respectively. The solid line is obtained from 
the model using Eq. 1161 with ro/(J — 1.0, and the dashed 
line is the liquid drop model with no heterogeneity. (Bot- 
tom) Critical barrier height, AF{n*)/kT, as a function of 
(1/lnS')^ Symbols are the same as above. 



IV. DISCUSSION 

One of the goals of this work was to develop an un- 
derstanding of how the presence of a small, nanoscale 
heterogeneity effects the condensation of a vapor and to 
quantify this by measuring both the nucleation rate and 
calculating the free energy barrier for the process. The 
obvious result is that the presence of the heterogeneity 
allows nucleation to occur faster, but when the attraction 
between the seed and the vapor is only marginally more 
attractive than the vapor-vapor interaction, we see the 
work of forming small clusters is essentially unchanged 
and the free energy is only lowered for the larger clus- 
ters. This appears to result from an interesting balance 
between the energetics of cluster formation and the en- 
tropic, or translational degrees of freedom, inherent in 
the system. There is a small energetic advantage to 
forming a monomer liquid particle that includes the seed, 
but there are many more ways of forming the monomer 
without the seed and these dominate the statistics in the 
density of states. Hence, there is no change in AF{n) 
for the small clusters. Larger clusters are rarer and the 
energetic advantage of including the seed in a cluster in- 
creases because it can have more neighbors. Eventually, 
the additional energy contribution dominates so the seed 
is connected to all the larger clusters and the barrier is 
lowered compared to homogeneous nucleation. 

We could completely decouple the heterogeneous and 
homogeneous processes by considering an order parame- 
ter that focuses just on n-sized clusters that include the 
seed. However, Fig. [5] suggests the height of the nucle- 
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ation barrier and the rate should remain unaffected be- 
cause the seed is always part of the cluster that goes 
over the barrier and it is part of this cluster long be- 
fore it reaches the critical size. As a result, the MFPT 
of the larger clusters will not change and neither should 
the probability of finding the larger clusters within the 
ensemble of runs. The work of forming the smaller clus- 
ters would be expected to increase, as the probability of 
finding monomers and dimers etc decreases relative to 
the homogeneous case and the probability of finding no 
liquid-like atoms connected to the seed increases. 

We also wished to make some quantitative assessment 
of the affect small heterogeneities may have on our abil- 
ity to measure homogeneous nucleation. The presence of 
a single seed in the current simulations changes the rate 
of nucleation by one or two orders of magnitude, but it 
is important to note that we have relatively high con- 
centration of heterogeneities in our system and we would 
expect, on the basis of the entropic arguments presented 
above, that the effect should decrease as the concentra- 
tion of heterogeneities decreases. 

Finally, the comparison of our simulation results with 
those of our model show that the simple capillarity base 
approach suffers the same degree of failure in heteroge- 
neous nucleation as it does in the homogeneous case. The 



models are able to make relatively good predictions con- 
cerning the size of the critical nucleus but not the height 
of the barrier. Nevertheless, the model does capture 
the general features of heterogeneous nucleation, demon- 
strating both a spinodal limit in the supersaturation for 
a fixed sized heterogeneity and a critical size of activa- 
tion, where the barrier to nucleation goes to zero for a 
completely wetting system above a certain size hetero- 
geneity. There are also many corrections that could be 
introduced to the model to account for the small size 
of the system, such as size dependent surface tensions 
or the disjoining pressure, which provides a thermody- 
namic description of the interaction of the interfaces in a 
very thin film, that will probably lead to some improve- 
ment in the agreement. However, the very concepts of 
surface wetting and contact angle become poorly defined 
when describing processes involving heterogeneities that 
are of the small size considered here and more molecular 
approaches may be needed to describe nanoscale hetero- 
geneous nucleation accurately. 
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